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Amendments to the Claims: 

Claim 1 . (Canceled) 

2. (Currently amended) The method of claim 3 4-, wherein the step of computing 
the global statistics further comprises the steps of estimating the global noise standard 
deviation a to generate the global statistics. 

3. (Currently Amended) Th e m e thod of claim 1, wh e rein the st e p of computing 
th e local statistics for each pix e l furth e r includes th e st e ps of: A method for reducing 
noise in a digital image formed from a plurality of pixels including a given pixel the 
method comprising the steps of: 

computing global statistics from the image: 
computing local statistics for the given pixel; 
configuring a local filter using the local and 
global statistics: 

filtering the given pixel using the local filter to reduce image noise: 

wherein the step of computing the local statistics for the given pixel 

further includes the steps of: 

selecting a window containing the given said pixel and a plurality of 
neighboring pixels; 

computing a the 2-D local variance of the given said pixel based on 
information related to the pixels in the window; 



3 of 23 



Patent Application No. 10/697,362 Docket No. SAM2.PAU.33 

Amdt. Dated February 14, 2007 

Reply to Office Action of November 20, 2006 

computing a plurality of t he 1-D local variances along multiple directions 
through the given said pixel in within the window; and 

detecting a the local edge direction by selecting one of the direction with 
the smallest 1-D local variance. 



4. (Currently amended) The method of claim 3 4-, wherein the step of computing 
the local statistics for the given pixel, e ach pix e l further includes the steps of: 

selecting a window containing the given said pixel and a plurality of 
neighboring pixels; 

computing the 2-D local variance <7q of the given said pixel based on 
information related to the pixels in the window; 

computing the 1-D local variances a] , a\ , or] , and <r 4 2 along the 

horizontal (Li), vertical (L 2 ), diagonal from upper left to lower right (L 3 ), and diagonal 
from upper right to lower left (L4) directions through the given said pixel, respectively, in 
within the window; and 

detecting the local edge direction by selecting on e of the direction 
dir e ction s with the smallest 1-D local variance. 



5. (Currently amended) The method of claim 3 4-, wherein the step of configuring 
the local filter for e ach pix e l using the local and global statistics further includes the steps 
of: 

selecting the detected local edge direction L as the direction of the local 

filter; 
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for the detected local edge direction L computing the 1 -D filter strength as 
a function of the square root of the local variance and the global noise standard deviation; 

computing the 2-D filter strength as a function of the local variance and 
the global noise standard deviation; and 

configuring the local filter for the detected local edge direction L based on 
the 1-D and 2-D filter strengths. 

6. (Currently amended) The method of claim 5, wherein the step of configuring 
the local filter for e ach pix e l using the local and global statistics further includes the steps 
of: 

selecting the detected local edge direction L* (k = 1, 2, 3, or 4) as the 
direction of the local filter; 

for the detected local edge direction L k computing the 1-D filter strength 
a k = min(2<T, max(3cr - a k ,0)) /(2<r) ; 

computing the 2-D filter strength a 0 = min(2<r, max(3<r -<7 0 ,0))/(2cr) ; 

and 

configuring the local filter f k for the detected local edge direction Lk 

according to the following conditions , wherein cu, is the filter strength along the local 
edge direction U fk = 0 (non-edge\ 1 (horizontal), 2 (verticals 3 (upper left to lower 
right, 4 (upper right to lower left)]: 

(i) if the detected direction is Li, then /, is configured as a 2-D local 
filter for horizontal direction, wherein: 
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•>> 9 



a 0 + 3a,(l-a 0 ) a 0 + 3(3-2a,)(l -ar 0 ) a 0 +3a,(l-a 0 ) 



a. 



(ii) if the detected direction is L 2 , then f 2 is configured as a 2-D local 
filter for vertical direction, wherein: 



h 9 



a 0 a 0 +3a 2 (l-a 0 ) a 0 
a 0 a 0 +3(3-2or 2 )(l-a 0 ) a 0 
a 0 a 0 +3or 2 (l-a 0 ) a 0 



(iii) if the detected direction is L 3 , then f 3 is configured as a 2-D local 
filter for the diagonal direction from upper left to lower right, wherein: 



3 9 



Qr 0 +3(3-2a 3 )(l-a 0 ) 



a 0 +3or 3 (l-a 0 ) 



; and 



(iv) if the detected direction is L 4 , then / 4 is configured as a 2-D local 
filter for the diagonal direction from upper right to lower left, wherein: 



f= l - 
U 9 



a* 



Qr 0 +3a 4 (l-a 0 )' 



a 0 + 3(3-2a 4 )(l-a 0 ) 



jx 0 +3or 4 (l-a 0 ) 



7. (Currently amended) The method of claim 2, wherein the steps of estimating 
the global noise standard deviation a further includes the steps of: 

dividing the input image into overlapping or non-overlapping blocks; 
computing a the mean and a the standard deviation for each block; 
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finding the smallest standard deviation d Q _ and its corresponding mean 
detecting block saturation due to noise; 

compensating for the smallest standard deviation^ to generate a 
compensated smallest standard deviation d 0 ; 

selecting the block standard deviations d n _ that are within a range of the 
compensated smallest standard deviation d 0 _; and 

averaging the selected block standard deviations d n to generate an 
estimate of the global noise standard deviation <j . 

8. (Currently amended) The method of claim 7, wherein the block size is 7 _x _ 7 
or 5 x 9 _ pixels. 

9. (Currently amended) The method of claim 7, wherein the steps of detecting 
saturation and compensating the smallest standard deviation further include the steps of 
d e t e rmining th e following : 

defining is an upper pixel value limit UL, a lower pixel value limit LL, and 
mid value M between UL and LL, 

wherein if the mean m 0 is less than the mid value rang e M 9 and the 

smallest standard deviation is greater than a the difference between the mean m 0 __ and the 
lower limit LL, then saturation has occurred at the lower limit LL, and the smallest 
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standard deviation d 0 is compensated by adding thereto a compensation term that is a 
function of the smallest standard deviation d 0 and said difference between the mean m 0 
and the lower limit LL, to generate the compensated smallest standard deviation d 0 ; 

else if the mid value range Mis less than the mean m Q , and the smallest 
standard deviation d 0 is greater than the difference between the upper limit UL and the 
mean m 0 , then saturation has occurred at the upper limit UL, and the smallest standard 
deviation d 0 is compensated by adding thereto a compensation term that is a function of 
the smallest standard deviation d Q and the said difference between the upper limit UL 
and the mean m 0 , to generate the compensated smallest standard deviation d 0 ; 

otherwise, no saturation has occurred, wherein d 0 =d 0 . 

10. (Currently amended) The method of claim 7, wherein the steps of detecting 
saturation and compensating the smallest standard deviation further include the steps of 
determining the following: 

where UL is an upper pixel value limit, LL is a lower pixel value limit, M 
is a mid pixel value, and UL<M< LL, if the mean m Q < M and the smallest standard 
deviation d 0 >m 0 -LL, then saturation has occurred at the lower limit LL, wherein d Q is 
compensated as d 0 = d Q + K • (d 0 - (m 0 - LL)) , such that AT is a compensation factor; 
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else if the mean m 0 >M and the smallest standard deviation 
d 0 >UL-m 0 , then saturation has occurred at the upper limit UL, wherein d 0 is 
compensated as d 0 = d 0 + K • (d 0 - (UL -m 0 )); 

otherwise, no saturation has occurred, wherein d 0 -d 0 . 

11. (Previously presented) The method of claim 10, wherein LL = 0,UL = 255, 
and M= 128. 

12. (Previously presented) The method of claim 7, wherein the step of selecting 
the block standard deviations further includes the steps of selecting the block standard 

deviation d n for averaging if | d n - d 0 |< max(<7 0 ,l) . 
Claim 13. (Canceled) 

14. (Currently amended) The system of claim 15 44, wherein the global statistics 
module estimates athe global noise standard deviation a to generate the global statistics. 

15. (Currently amended) Th e system of claim 13, A noise reduction system for 
reducing noise in a digital image comprising pixels, the system comprising: 

a global statistics module that computes global statistics from the image; 
a local statistics module that computes local statistics for each of a 
plurality of image pixels including a given pixel: 
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a filter configuration module that uses the local and global statistics for 
the given pixel to configure a local filter for filtering the given pixel: 

the local filter as configured by the filter configuration module, being 

adapted for filtering the given pixel to reduce image noise , wherein the local statistics 
module computes the local statistics for the given eaeh pixel by: 

selecting a window containing the given said pixel and a plurality of 
neighboring pixels; 

computing a the 2-D local variance of said pixel based on information 
related to the pixels in the window; 

computing a plurality of the 1-D local variances along multiple directions 
each defined by a pair of the pixels through said pix e l in within the window: and 

detecting a the local edge direction for the given pixel by selecting on e of 
the direction with the smallest 1-D local variance. 



16. (Currently amended) The system of claim 15 44, wherein the local statistics 
module computes the local statistics for each pixel by: 

selecting a window containing said pixel and a plurality of neighboring 

pixels; 

computing the 2-D local variance a] of the given sakl pixel based on 
information related to the pixels in the window; 

computing the 1-D local variances a] , o\ , a\ , and a\ along the 
horizontal (Li), vertical (L 2 ), diagonal from upper left to lower right (L3), and diagonal 
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from upper right to lower left (L 4 ) directions through the given said pixel, respectively, in 
within the window; and 

detecting the local edge direction by selecting efte-ef the direction 
dir e ctions with the smallest 1-D local variance. 

17. (Currently amended) The system of claim 15 43, wherein the filter 
configuration module configures the local filter for each pixel using the local and global 
statistics by: 

selecting a the detected local edge direction L as the direction of the local 

filter; 

for the detected local edge direction L computing a the 1-D filter strength 
as a function of the square root of the local variance and the global noise standard 
deviation; 

computing a the 2-D filter strength as a function of a the local variance 
and a the global noise standard deviation; and 

configuring the local filter for the detected local edge direction L based on 
the 1-D and 2-D filter strengths. 

18. (Currently amended) The system of claim 17, wherein the filter configuration 
module configures configured the local filter for each pixel using the local and global 
statistics, wherein c k is the filter strength and at is the global noise standard deviation for 
edge directions U \ k = 0 fnon-edee\ 1 (horizontal!. 2 (vertical^). 3 (upper left to lower 
right). 4 (upper right to lower left)] bv: 
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selecting the detected local edge direction L* (k = 1, 2, 3, or 4) as the 
direction of the local filter; 

for the detected local edge direction I* a 

computing the 1 -D filter strength a k = min(2cr, max(3<r - a k ,0)) /(2a) A 

and [[;]] computing the 2-D filter strength 



according to the following conditions: 

(i) if the detected direction is Li, then f } is configured as a 2-D local 
filter for horizontal direction, wherein: 



(ii) if the detected direction is L 2 , then f 2 is configured as a 2-D local 
filter for vertical direction, wherein: 



(iii) if the detected direction is L 3 , then f 3 is configured as a 2-D local 
filter for the diagonal direction from upper left to lower right, wherein: 



a 0 = min(2<r,max(3o' -<r 0 ,0))/(2<7) ; and 



configuring the local filter f k for the detected local edge direction Lk 





a 0 a 0 +3<* 2 (l-a 0 ) a 0 
a 0 a 0 + 3(3-2tf 2 )(l-tf 0 ) a 0 
a Q a 0 +3a 2 (\-a Q ) a 0 



tf 0 + 3or 3 (l-a 0 ) 




a 0 + 3« 3 (l-a 0 ) 




; and 
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(iv) if the detected direction is L4, then f 4 is configured as a 2-D local 
filter for the diagonal direction from upper right to lower left, wherein: 



u 9 



a 0 +3a 4 (l-a 0 j 



a 0 +3(3-2a 4 )(l-or 0 ) 



a 0 +3a 4 (l-a 0 ) 



19. (Currently amended) The system of claim 14, wherein the input image 
comprises a plurality of ov e rlapping or non ov e rlapping blocks, and wherein the global 
statistics module further comprises: 

a mean and standard deviation module that computes the mean and the 
standard deviation for each block ; 

a minimum finder module that finds the smallest standard deviation d 0 and 
its corresponding mean m 0 ; 

a saturation detector that detects block saturation due to noise; 

a saturation compensator that compensates for the smallest standard 
deviation^ to generate a compensated smallest standard deviation 3 0 ; and 

a selective averaging module that selects the block standard deviations 
d n that are within a range of the compensated smallest standard deviation d 0 , and 
averages the selected block standard deviations d„ to generate an estimate of the global 
noise standard deviation a . 
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20. (Currently amended) The system of claim 19, wherein the block size 
is7_x_7 or 5x9_pixels. 

21. (Currently amended) The system of claim 19, wherein: 

an upper pixel value limit is denoted UL, a lower pixel value limit is 
denoted LL, and a mid value M is between UL and LL 9 

wherein the saturation detector determines if the mean m 0 is less than the 
mid range M, and the smallest standard deviation is greater than a the difference between 
the mean m 0 and the lower limit II, indicating that saturation has occurred at the lower 
limit LL, and if so, the saturation compensator compensates the smallest standard 
deviation d 0 [[is]] by adding thereto a compensation term that is a function of the 

smallest standard deviation d 0 and said difference between the mean m 0 and the lower 

limit LL, to generate the compensated smallest standard deviation d 0 ; 

else if the saturation detector determines that the mid value range Mis less 
than the mean m 0 , and the smallest standard deviation d 0 is greater than a the difference 

between the upper limit UL and the mean m 0 , indicating saturation has occurred at the 
upper limit UL, the saturation compensator compensates the smallest standard deviation 
d 0 by adding thereto a compensation term that is a function of the smallest standard 

deviation d 0 and said difference between the upper limit UL and the mean m Q , to 
generate the compensated smallest standard deviation d 0 ; 

otherwise, no saturation has occurred, wherein d 0 =d Q . 
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22. (Previously presented) The system of claim 21 , wherein: 

if the mean m 0 <M and the smallest standard deviation d Q >m 0 -LL , 

indicating saturation has occurred at the lower limit LL, then d 0 is compensated as 
d 0 =d 0 + K -(d 0 - (m 0 - LL)) , such that K is a compensation factor; 

else if the mean m 0 >M and the smallest standard deviation 
d 0 >UL-m 0 , indicating saturation has occurred at the upper limit UL, then d Q is 
compensated as d Q = d 0 + K • (d Q - {UL - m Q )) ; 

otherwise, no saturation has occurred, wherein d 0 =d 0 . 

23. (Previously presented) The system of claim 22, wherein LL = 0 9 UL = 255, 
andM= 128. 

24. (Previously presented) The system of claim 19, wherein block standard 
deviations d n are selected for averaging if | d n - d 0 \< max(d 0 ,l) . 
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